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strongly correlated electron system as a function of temperature and interac- 
tion strength. Within a dynamical mean field theory of the Hubbard model 
that becomes exact in the limit of large lattice coordination, we predict an 
anomalous enhancement of spectral weight as a function of temperature in 
the correlated metallic state and report on experimental measurements which 
agree with this prediction in V 203 - We argue that the optical conductivity 
anomalies in the metal are connected to the proximity to a crossover region 
in the phase diagram of the model. 
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The interest in the distribution of spectral weight in the optical conductivity of correlated 
electron systems has been revived by the improvement in the quality of the experimental 
data [|^Q. The traditional theoretical methods used in the strong correlation problem, 
have only been partially successful in describing the interesting behavior of the optical 
response which takes place in the strong correlation regime. This is most notable when the 
temperature dependence of the spectral weight is considered. We present here new optical 
conductivity data on metallic V2O3 and analyze the problem theoretically it using a mean 
held theory which is exact in the limit of large number of spatial dimensions P]. In this limit, 
lattice models can be mapped onto an equivalent impurity model subject to a selfconsistency 
condition for the conduction electron bath |Q,|^ . This technique has already given some new 
insights into the Mott transition, which is realized in the solution of the Hubbard model 
0 -^. The new experimental measurements on V2O3 test a critical prediction of the theory. 

The goal of this paper is to demonstrate two points: i) that the present dynamical mean 
held approach is a useful tool which allows us to make qualitative predictions such as an 
anomalous enhancement of spectral weight as a function of temperature. This behavior 
of the spectral weight is conhrmed by experiments in V2O3 reported here; and ii) that 
the Hubbard model on a frustrated lattice, in the limit of large lattice coordination, with 
parameters extracted from the optical data, accounts for many other experimental features 
of the V2O3 system, such as its phase diagram and large specihc heat capacity. We also 
present results for the various quantities which are useful in the analysis of the optical 
conductivity of the Hubbard model , treated in an approximation which is exact in a well 
dehned limit. 

The experiments were carried out on single crystals of V2O3 which were grown by the 
skull-melting process, carefully annealed, polished and mounted on a temperature-controlled 
stage in an optical cryostat |]I|,pT|. The specular reflectivity was measured using a Fourier 
transform spectrometer over a frequency range from 6meV to 3 . 7 eH. The data were com- 
bined with dc conductivity and higher energy reflectivity measurements, and a Kramers- 
Kronig transformation was used to determine the conductivity. 
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The theoretical model is the Hubbard hamiltonian, 



H = - Y. + YU (1) 

<i,j> i ^ ^ 

where summation over repeated spin indices is assumed. 

In the limit of inhnite dimensions the coordination number of the lattice (i.e. the number 
of neighbors g) gets large and the hopping is scaled as tjj — H]- The lattice model is 
mapped onto an equivalent impurity problem supplemented by a selfconsistency condition. 
The derivation of this result can be found elsewhere 0,^, so we only present the hnal 
expressions. The resulting local effective action reads, 



Siocai = - [ dr ( dr'4(r)G'o^(r - r')c^(r') + U f dT{n^{r) - ^)(n|(r) - (2) 

Jo Jo Jo 2 2 



where cj., correspond to a particular site. Requiring that Giocaiij^) = ^kG{k,uj) the 
selfconsistency condition becomes Gq = a; + /r — t'^Giocaii^), where we have assumed a 
semi-circular bare density of states p°{e) = {2/7iD)^Jl — {e/Dy, with t = y, which can be 
realized in a Bethe lattice and also on a fully connected fully frustrated version of the model 
m- We consider the symmetric case /i = 0. We use an exact diagonalization algorithm 
(ED) and an extension of the second order perturbative (20PT) calculation to solve the 
associated impurity problem [p]] . 

Early work on the optical conductivity in the Hubbard model was carried out by Pruschke 
et al. using quantum Monte Carlo 0, a technique that is restricted to high temperatures. 
The 20PT and ED techniques allow us to study the low temperature regime near the Mott 
transition relevant to the experiments addressed here. 

The optical conductivity is defined as cr{uj) = — ^Im([j, j]). In the limit of d = cxo it can 
be expressed in terms of the one particle spectrum |7|JT3[|, 

, , 1 2eH^a^ dio' 



a[oo = 



/ OO PCO rj( y 

dep°{e) / -Ayj)Ayj + u;){nf{J)-nf{J + u;)) 

-OO J —OO 2i7T 



(3) 



oj vTi^ 

with A^ui) = — 2Im[Gfc(o;)] , e the electron charge, z/ the volume of the unit cell, and a the 
lattice constant. 
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At T = 0, <j{u) can be parametrized by <j{uj) = + aregi^) where the coefficient 

in front of the 5-fnnction is the Drnde weight and ujp is the renormalized plasma freqnency 



In the presence of disorder S{u!) is replaced by a lorentzian of width F. The expectation 



valne of the kinetic energy {K) is related to the condnctivity by the snm rnle 



roc 

/ a{u)du; = = 
Jo 2dh V 



UJ p 

Att 



( 4 ) 



In the limit d — >■ cx) the Drnde weight is obtained in terms of the qnasiparticle weight Z, 



4tt 



4-7rt^e^Q^ 

FT 



Zp“(0). 



The solntion of the mean held eqnations shows that at low temperatnres the model 
has a metal insnlator transition (Mott-Hnbbard transition) at an intermediate valne of the 
interaction Uc ~ 3D p-^]. The metallic and insnlator solntions at low T are very different 
and we schematically represent them in Fig. |l[ The optical condnctivity response, in a hrst 
approximation, can be nnderstood from transitions between the states which appear in the 
the DOS. We schematically sketch the optical response in the lower part of Fig. 

We will discnss the resnlts for the model in regard of different experimental data on 
the V 2 O 3 system. Vanadinm oxide has three t 2 g orbitals per V atom which are hlled with 
two electrons. Two electrons (one per V) are engaged in a strong cation-cation bond, 
leaving the remaining two in a twofold degenerate Cg band |^. LDA calculations give a 
bandwidth of ~ 0.5eV [|T^]. The Hubbard model ignores the degeneracy of the band which 



is crucial in understanding the magnetic structure [^, but captures the interplay of the 
electron-electron interactions and the kinetic energy. This delicate interplay of itinerancy 
and localization is responsible for many of the anomalous properties of this compound, and 
are correctly predicted by this simplihed model. 

Experimentally one can vary the ratio ^ by introducing V vacancies. The parameters 
are extracted from the experimental data on V 2 -yO^ by comparing the measurments with 
the schematic spectra displayed in Fig. |^, and are summarized in table It is not surprising 
that U and D are different in the metal and the insulator, since the lattice parameter and 
the screening length change rapidly across the phase boundary. 
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Our calculations have been performed on the model with nearest neighbor hopping 
and next to nearest neighbor hopping ^ on a Behte lattice. The n.n.n. hopping introduces 
magnetic frustration which is essential if we want to describe V 2 O 3 with a one band model. 
The condition t\ keeps the bare density of states invariant. For ^ 2/^1 = 0 we 

recover the original hamiltonian, and ^ 2/^1 = 1 gives the PM solution. The same model on 
the hypercubic lattice would give very similar results. In Fig. we display the phase diagram 
for = y'^ 1/3. It has the same topology as the experimental one Ill0| , [l7|, [l^ . Frustration 
lowers the T^^ei well below the second order Tmit point [§]. Using the parameters of table 
|, Tmit ~ 240iF, which is only within less than a factor of 2 from the experimental result. 
The dotted line indicates a crossover separating a good metal at low T and a semiconductor 
at higher T. Between these states pdc has an anomalous rapid increase with T. The reason 
for this feature can be traced to the disappearance of the coherent central quasiparticle peak 
in the DOS. We find the behavior of pdc to be in good agreement with the experimental 



results of Me Whan et al. (inset Fig. [^]. Another crossover is indicated by a shaded 
area, it separates a semiconducting region with a gap A comparable with T from a good 
insulator where T << A, consequently the crossover temperature increases linearly with U 
and the crossover width becomes broader with increasing T. This crossover is characterized 
by a sudden increase in pdc as function of U at a fix T (inset Fig. ||). This crossover was 
experimentally observed in V 2 O 3 by Kuwamoto et al. 0 - 



The experimental optical spectrum of the insulator is reproduced in Fig. || (dotted lines). 
It is characterized by an excitation gap at low energies, followed by an incoherent feature 
that corresponds to charge excitations of mainly Vanadium character |I|. These data are to 
be compared with the model results of Fig. The size of the gaps A (shown for various 
degrees of magnetic frustration in the inset of Fig. |^), and the overall shape of the spectrum 
are found to be in good agreement with the experimental results []^. Another quantity that 

J2 

can be compared to the experiment is the integrated spectral weight which is related to 
{K) by the sum rule (^). Setting the lattice constant a ~ 3 A the average V — V distance, we 
hnd our results also in good agreement with the experiment (left inset of Fig. ^). Notice that 
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the dependence of {K) and A on the degree of magnetic frustration is not very strong, which 
is consistent with the view that magnetism in this system is not the driving mechanism of 
the metal insulator transition. Notice also that the experimental points are closer to the 
intermediate frustration curves than to the unfrustrated ones. 

We now discuss the new data in the metallic phase. The experimental data for pure 
samples that become insulating at ~ 150W were obtained for T = 170K and T = 300W, 
and are presented in Fig. || (full lines). Both spectra are made up of broad absorption at 
higher frequencies and some phonon lines in the far infrared. They appear to be rather 
featureless, however, upon considering their difference (in which the phonons are approxi- 
mately eliminated) distinct features are observed. As T is lowered, there is an enhancement 
of the spectrum at intermediate frequencies of order 0.5eV; and more notably, a sharp low 
frequency feature emerges that extends from 0 to O.lSel^. Moreover, these enhancements 

J2 

result in an anomalous change of the total spectral weight with T. We argue below, that 
these observations can be accounted by the Hubbard model treated in mean held theory. 

In Fig. ^ we show the calculated optical spectra obtained from 20PT for two different 
values of T. The interaction is set to U = 2. ID that places the system in the correlated 
metallic state. It is clear that, at least, the qualitative aspect of the physics is already 
captured. As the temperature is lowered, we observe the enhancement of the incoherent 
structures at intermediate frequencies of the order ^ to U and the rapid emergence of a 
feature at the lower end of the spectrum. Setting D pe OAeV we hnd these results consistent 
with the experimental data on V 2 O 3 . The two emerging features can be interpreted from 
the qualitative picture that was discussed before which is relevant for low T. 

An interesting prediction of the model is the anomalous increase of the integrated spectral 
weight ^ as T is decreased. This is due to the rather strong T dependence of the kinetic 
energy {K) oc ^ in the region near the crossover indicated by a dotted line in the phase 
diagram. This also corresponds with the behavior of the low frequency features in the 
spectra. From the model calculations we expect an enhancement of the spectral weight as 
T is decreased. With the chosen parameters this should occur at a scale T^oh ~ 0.05D ~ 
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2A0K which correlates well with the experimental data. Tcoh has the physical meaning of 

the temperature below which the Fermi liquid description applies 0, as the quasiparticle 

2 

resonance emerged in the density of states. We hnd ^ ~ which is somewhat lower 

than the experimental result. This could be due to the contribution from tails of bands at 
higher energies that are not included in our model, or it may indicate that the bands near 
the Fermi level are degenerate. 

Experiment show that the slope of the linear term in the specihc heat 7 in the metallic 
phase is unusually large. For 0.08 Ti substitution 7 40 ^^^^2 , while for a pressure of 

2bKbar in the pure compound 7 ~ with V dehciency in a range of y = 0.013 

to 0.033 the value is 7 ~ 47 ^^'^2 ||T^ . In our model 7 is simply related to the weight in the 
Drude peak in the optical conductivity and to the quasiparticle residue Z , 7 = ^3 . 

The values of U = 2.1D and D 0.4eE extracted from the optical data correspond to a 
quasiparticle residue Z k, 0.3, and results in 7 ‘^^ rrwiK^ which is close to the experimental 
hndings. Thus, it turns out that the mean held theory of the Mott transition explains in a 
natural and qualitative manner, the experimentally observed optical conductivity spectrum, 
the anomalously large values of the slope of the specihc heat 7, and the dc conductivity 
in the metallic phase as consequence of the appearance of a single small energy scale, the 
renormalized Fermi energy e*p. 

To conclude, we have illustrated how the mean held theory, that becomes exact in the 
limit of large dimensions, can be used to study the physics of systems where the interactions 
play a major role. We presented new data on the metallic phase of V2O3 revealing an 
anomalous enhancement of the optical spectral weight as the temperature is lowered in 
agreement with the surprising prediction of the mean held theory. We also contrasted several 
experimental results to the solution of the model Hamiltonian to argue that, in the mean 
held approximation the Hubbard model gives a consistent semiquant it at ive picture of this 
strongly correlated electron system. A great challenge is to extend the mean held approach 
to incorporate more realistic band structure density of states with orbital degeneracies and 
more complicated unit cells. This extensions will require more efficient tools for solving 
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larger impurity models, and would allow for a more quantitative description of the physical 
properties of transition metal oxide systems. 

We acknowledge valuable discussions with V. Dobrosavljevic, G. Moeller A. Ruckenstein 
Q. Si and C.M. Varma. This work was supported by the NSF grant DMR 92-24000 




REFERENCES 



* Present Address: LPS, Ecole Normale Superieure, 24 rue Lhomond, 75231 Paris Cedex, 
05 France. 

[1] G. A. Thomas et al., Phys. Rev. Lett. 73, 1529 (1994); G. A. Thomas et ai., Journ. 
Low Temp. Phys. 95 , 33 (1994). 

[2] B. Bucher, Z. Schlesinger, P. G. Ganheld and Z. Fisk, Phys. Rev. Lett. 72 , 522 (1994). 

[3] W. Metzner and D. Vollhardt, Phys. Rev. Lett. 62, 324 (1989). 

[4] A. Georges and G. Kotliar, Phys. Rev. B 45, 6479 (1992). 

[5] A. Georges, G. Kotliar and Q. Si, hit. J. Mod. Phys. B 6, 705 (1992). 

[6] M. Jarrell, Phys. Rev. Lett. 69 , 168 (1992); M. J. Rozenberg, X. Y. Zhang and G. 

Kotliar, Phys. Rev. Lett. 69 , 1236 (1992); A. Georges and W. Krauth, Phys. Rev. Lett. 
69 , 1240 (1992). 

[7] Th. Pruschke, D. L. Gox and M. Jarrell, Phys. Rev. B 47, 3553 (1993). 

[8] A. Georges and W. Krauth, Phys. Rev. B 48, 7167 (1993). 

[9] M. J. Rozenberg, G. Kotliar and X. Y. Zhang, Phys. Rev. B 49 , 10181 (1994). 

[10] S. A. Garter, T. F. Rosenbaum, J. M. Honig, and J. Spalek, Phys. Rev. Lett. 67 , 3440 
(1991). 

[11] X. Y. Zhang, M. J. Rozenberg and G. Kotliar, Phys. Rev. Lett. 70 , 1666 (1993). 

[12] M. Gaffarel and W. Krauth, Phys. Rev. Lett. 72 , 1545 (1994); Q. Si, M. Rozenberg, 

G. Kotliar and A. Ruckenstein, Phys. Rev. Lett. 72 , 2761 (1994); M. J. Rozenberg, G. 
Moeller, and G. Kotliar, Mod. Phys. Lett. B 8, 535 (1994). 

[13] Anil Khurana, Phys. Rev. Lett. 64, 1990 (1990). 

[14] W. Kohn, Phys. Rev. 133, A171 (1964). 



9 




[15] C. Castellani, C. R. Natoli, and J. Ranninger, Phys. Rev. B 18 , 4945 (1978). 

[16] L. F. Mattheiss, J. Phys., Condens. Matter 6, 6477 (1994). 

[17] H. Kuwamoto, J. M. Honig, and J. Appel, Phys. Rev. B 22 , 2626 (1980). 

[18] D. Me Whan et ah, Phys. Rev. B 7, 1920 (1973). 

[19] D. Me Whan et ah, Phys. Rev. Lett. 27 , 941 (1971); D. Me Whan et ah, Phys. Rev. B 
7 , 3079 (1973); S. A. Carter et ah, Phys. Rev. B 48 , 16841 (1993). 



10 




FIGURES 

FIG. 1. Schematic DOS ( 1/2 filling) and their corresponding optical spectra for the metallic 
and insulator solutions. The approximate width of the incoherent peaks in the DOS is 2D and 
that of the central peak in the metal is ZD = e*p. 

FIG. 2 . Approximate phase diagram for the model with n.n. and n.n.n. hopping 
(G/^i) = \/l/3. The 1st order PM metal-insulator transition ends at the critical point Tmit 
(square). The dotted line and the shaded region describe two crossovers as discussed in the text. 
The full circles indicate the position of the optical spectra. A: insulator (y = 0), B\ insulator 
(y = .013), C\ metal (y = 0, 170AT), D: metal (y = 0, 300AT). Note that for comparison with ex- 
perimental results increasing U/D is associated with decreasing pressure Left inset: Pdc{T) 

for U/D = 2.1, 2.3, 2.5 (bottom to top). The maxima of Pdc{T) defines the dotted line. Right inset: 
Pdc(U) for T = 0.06D (full) and T = 0.15D (dotted). 

FIG. 3. The experimental cr{u!) of metallic V 2 O 3 (full lines) at T = 170AT (upper) and T = 300iF 
(lower). The inset contains the difference of the two spectra Aa{u>) = anoKi^^) ~ o' 300 i<'(‘^)- 
Diamonds indicate the measured dc conductivity adc- Dotted lines indicate cr{u>) of insulating 
V 2 -y 03 with y = .013 at lOAT (upper) and y = 0 at 70iF (lower). 

FIG. 4. a{u>) for the metallic solution (full lines) at U = 2. ID and T = 0.05D (upper) and 
0.083D (lower). A small F = 0.3 and 0.5D was included to mimic a finite amount of disorder. 
Dotted lines indicate the insulating solution results at U = 4D and T = 0 from ED (thin-dotted) 
and 20PT (bold-dotted). The peaks in ED result from the finite size of the clusters that can be 
considered. Left inset: {K) versus U for the AEM (bold-dotted) and PM insulators (thin), PM 
metal (bold) and partially frustrated model (thin-dotted). Right inset: A versus U for the AFM 
(dotted), partially frustrated (thin) and PM (bold) insulators. A is twice the energy of the lowest 
pole from the ED Green function. The data are for NgUes 00 from clusters of 3, 5 and 7 sites 
assuming l/NgUes scaling behavior. Black squares show the insulator experimental results. 
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TABLES 



Phase 


D [eV] 


U [eV] 


Parameter 

A [eV] 


u;p/47r [eV/flcm] 


Insulator (y=0) 


0.33 ± 0.05 


1.3 ±0.05 


0.64 ± 0.05 


170 ± 20 


Insulator (y=.013) 


0.46 ± 0.05 


0.98 ± 0.05 


0.08 ± 0.05 


800 ± 50 


Metal (170K) 


0.4 ±0.1 


0.8 ±0.1 


- 


1700 ± 300 



TABLE I. Experimental parameters for the model. 
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